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Context. Evolution of spiral galaxies in the Virgo cluster. 

Aims. We determined the parameters of the interaction between the interstellar medium of NGC 4522 and the intracluster 
medium of the Virgo cluster. 

Methods. A dynamical model including ram pressure stripping is applied to the strongly Hi deficient Virgo spiral galaxy 
NGC 4522. A carefully chosen model snapshot is compared with existing VLA Hi observations. 

Results. The model successfully reproduces the large-scale gas distribution and the velocity field. However it fails to reproduce 
the large observed Hi linewidths in the extraplanar component, for which we give possible explanations. In a second step, we 
Q ' solve the induction equation on the velocity fields of the dynamical model and calculate the large scale magnetic field. Assuming 
^ I a Gaussian distribution of relativistic electrons we obtain the distribution of polarized radio continuum emission which is also 
■ compared with our VLA observations at 6 cm. The observed maximum of the polarized radio continuum emission is successfully 
reproduced. Our model suggests that the ram pressure maximum occurred only ~ 50 Myr ago. 

Conclusions. Since NGC 4522 is located far away from the cluster center (~ 1 Mpc) where the intracluster medium density is too 
T ""j ' low to cause the observed stripping if the intracluster medium is static and smooth, two scenarios are envisaged: (i) the galaxy 
I^nJ ' moves very rapidly within the intracluster medium and is not even bound to the cluster; in this case the galaxy has just passed 
the region of highest intracluster medium density; (ii) the intracluster medium is not static but moving due to the infall of the 
M49 group of galaxies. In this case the galaxy has just passed the region of highest intracluster medium velocity. This study 
shows the strength of combining high resolution Hi and polarized radio continuum emission with detailed numerical modeling 
of the evolution of the gas and the large-scale magnetic field. 

Key words. Galaxies: individual: NGC 4522 - Galaxies: interactions - Galaxies: ISM - Galaxies: kinematics and dynamics 



Oh 



1. Introduction 

The spiral galaxy NGC 4522 is one of the best examples for 
ongoing ram pressure stripping due to the galaxy's rapid 
motion within the hot and tenuous intracluster gas (ICM) 
of the Virgo cluster. Hi and Ha observations (Kenney et al. 
2004, Kenney & Koopmann 1999) showed a heavily trun- 
cated gas disk at a radius of 3 kpc, which is ~ 40% of the 
optical radius, and a significant amount of extraplanar gas 
to the west of the galactic disk. The one-sided extrapla- 
nar atomic gas distribution shows high column densities, 
comparable to those of the adjectant galactic disk. Such 
high extraplanar gas column densities are unusual among 
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Virgo spiral galaxies (Cayatte et al. 1990), suggesting that 
NGC 4522 is in a short-lived phase of evolution. 

Since the stellar disk is symmetric and undisturbed 
(Kenney & Koopmann 1999), a tidal interaction is ex- 
cluded as the origin of the peculiar gas distribution of 
NGC 4522. A scenario where NGC 4522 is moving roughly 
to the east and experiences strong ram pressure can ac- 
count for the truncated gas disk and the western extra- 
planar gas. A scenario where the ram pressure peak oc- 
curred a few 100 Myr ago (stripping by the intracluster 
medium around M87; VoUmer et al. 2000) is excluded, 
because the Hi column density distribution and velocity 
field are inconsistent with fall back of stripped gas which 
should occur a few 100 Myr after peak ram pressure, which 
is needed to create extraplanar gas at these late times. 
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Further evidence for the peak ram pressure scenario comes 
from polarized radio continuum observations (VoUmer et 
ah 2004). The 6 cm polarized emission is located at the 
eastern edge of the galactic disk, opposite to the western 
extraplanar gas. This ridge of polarized radio continuum 
emission is most likely due to ram pressure compression 
of the interstellar medium (ISM) and its magnetic field. 
In addition, the degree of polarization decreases from the 
east to the west and the flattest spectral index between 
20 cm and 6 cm coincides with the peak of the 6 cm po- 
larized emission. These findings are also consistent with a 
scenario where ram pressure is close to its maximum. 

This scenario has one important caveat: NGC 4522 is 
located at a projected distance of '--^ 1 Mpc from the cen- 
ter of the Virgo cluster (M87). Assuming a static smooth 
ICM and standard values for the ICM density and the 
galaxy velocity, the ram pressure at that location seems 
to be too low by an order of magnitude to produce the 
observed truncation of the gas disk. Therefore, Kenney et 
al. (2004) and VoUmer et al. (2004) claim that either (i) 
NGC 4522 has a velocity of 4000 kms"^ with respect 
to the Virgo cluster mean or (ii) that the ICM density is 
locally enhanced or (iii) that the intracluster medium is 
moving or (iv) a combination of (ii) and (iii). Based on ex- 
isting X-ray observations Kenney et al. (2004) claim that 
the latter scenario is more likely, given that NGC 4522 is 
located between M87 and the bright elliptical galaxy M49 
which is the center of a galaxy subcluster falling into the 
Virgo cluster from behind (Irwin & Sarazin et al. 1996 and 
Biller et al. 2004). In addition, the X-ray hotspot detected 
by Shibata et al. (2001) represents the sign of an interac- 
tion between the Virgo intracluster medium and that of 
the M49 subcluster. The intracluster medium of the M49 
subcluster has thus a velocity whose radial component is 
opposite to the radial velocity of NGC 4522. This moving 
intracluster medium together with a modest local density 
enhancement of the intracluster medium by a factor of 2-4 
can account for the observed stripping radius. 

In this article we present a ram pressure stripping 
model for NGC 4522 in order to assess if a simulation 
using a time dependent ram pressure can account for the 
observed gas distribution, velocity field, and polarized ra- 
dio continuum emission distribution. Our aim is to con- 
strain the parameters of the event. Since the observations 
of Kenney et al. (2004) and VoUmer et al. (2004) showed 
that ram pressure is ongoing, we only consider model snap- 
shots close to peak ram pressure. 

In Sect.[5]we describe the dynamical model. The choice 
of the best fit model is justified and its time evolution and 
the final snapshot for the comparison with observations is 
presented in Sect. 13 The comparison with VLA Hi data is 
done using moment maps and position-velocity diagrams 
(Sect. 2)). We solved the induction equation for the veloc- 
ity field of our simulation to obtain the evolution of the 
large scale magnetic field and the distribution of the po- 
larized radio continuum emission which is compared with 
observations in Sect. |S1 We discuss our results in Sect. 
and give our conclusions in Sect. [3 



Table 1. Total mass, number of particles N, particle mass 
M, and smoothing length I for the different galactic com- 
ponents. 



component 


Mtot (Mq) 


N 


M (Mo) / (pc) 


halo 


4.8 10'" 


32768 


1.48 10" 1200 


bulge 


1.710^ 


16384 


1.010^ 180 


disk 


8.310® 


32768 


2.5 lO'^ 240 



2. Dynamical model 

We use the N-body code described in VoUmer et al. (2001) 
which consists of two components: a non-coUisional com- 
ponent that simulates the stellar bulge/disk and the dark 
halo, and a coUisional component that simulates the ISM. 

The non-coUisional component consists of 49 125 par- 
ticles, which simulate the galactic halo, bulge, and disk. 
The characteristics of the different galactic components 
are shown in Table Q The resulting rotation velocity is 
^100 kms^^ and the rotation curve becomes flat at a ra- 
dius of about 4 kpc. 

We have adopted a model where the ISM is simulated 
as a coUisional component, i.e. as discrete particles which 
possess a mass and a radius and which can have inelas- 
tic collisions (sticky particles). Since the ISM is a turbu- 
lent and fractal medium (see e.g. Elmegreen & Falgarone 
1996), it is neither continuous nor discrete. The volume 
filling factor of the warm and cold phases is smaller than 
one. The warm neutral and ionized gas fill about 30 — 50% 
of the volume, whereas cold neutral gas has a volume fill- 
ing factor smaller than 10% (Boulares & Cox 1990). It 
is not clear how this fraction changes, when an external 
pressure is applied. In contrast to smoothed particles hy- 
drodynamics (SPH), which is a quasi continuous approach 
and where the particles cannot penetrate each other, our 
approach allows a finite penetration length, which is given 
by the mass-radius relation of the particles. Both methods 
have their advantages and their limits. The advantage of 
our approach is that ram pressure can be included eas- 
ily as an additional acceleration on particles that are not 
protected by other particles (see VoUmer et al. 2001). In 
this way we avoid the problem of treating the huge density 
contrast between the ICM (n ^ 10~* cm~^) and the ISM 
{n > 1 cm~'^) of the galaxy. 

The 20 000 particles of the coUisional component repre- 
sent gas cloud complexes which are evolving in the gravita- 
tional potential of the galaxy. The total assumed gas mass 
is Mg°* = 1.3 10^ Mq, which corresponds to the total neu- 
tral gas mass before stripping, i.e. to an Hi deficiency of 0, 
which is defined as the logarithm of the ratio between the 
Hi content of a field galaxy of same morphological type 
and diameter and the observed Hi mass. To each parti- 
cle a radius is attributed depending on its mass. During 
the disk evolution the particles can have inelastic colli- 
sions, the outcome of which (coalescence, mass exchange. 
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or fragmentation) is simplified following Wiegel (1994). 
This results in an effective gas viscosity in tlie disk. 

As tlie galaxy moves tlirough the ICM, its clouds are 
accelerated by ram pressure. Within the galaxy's inertial 
system its clouds are exposed to a wind coming from a 
direction opposite to that of the galaxy's motion through 
the ICM. The temporal ram pressure profile has the form 
of a Lorentzian, which is realistic for galaxies on highly 
eccentric orbits within the Virgo cluster (Vollmer et al. 
2001). The effect of ram pressure on the clouds is sim- 
ulated by an additional force on the clouds in the wind 
direction. Only clouds which are not protected by other 
clouds against the wind are affected. 

The particle trajectories are integrated using an adap- 
tive timestep for each particle. This method is described 
in Springel et al. (2001). The following criterion for an 
individual timestep is applied: 



ai 

where is the acceleration of the particle i. The mini- 
mum value of ti defines the global timestep used for the 
Burlisch-Stoer integrator that integrates the coUisional 
component. 

3. Search for the best fit model 

In this section we will constrain the parameters of the ram 
pressure stripping event which are (i) the peak ram pres- 
sure, (ii) the temporal ram pressure profile, (iii) time since 
peak ram pressure, (iv) the inclination angle between the 
galaxy's disk and the intracluster medium wind direction, 
and (v) the azimuthal viewing angle for the observed incli- 
nation and position angles. These parameters are related 
to the observed quantities which are (1) the position an- 
gle, (2) the inclination angle of the galactic disk, (3) the 
line-of-sight velocity of the galaxy with respect to the clus- 
ter mean, and (4) the projected ICM wind direction. The 
position angle and inclination of NGC 4522 define a plane 
in three dimensional space. The model galaxy can then be 
rotated within this plane by the azimuthal viewing angle 
(see below). The three dimensional model wind direction, 
the line-of-sight velocity of the galaxy, and the projected 
ICM wind direction are thus functions of the azimuthal 
viewing angle. 

In the case of a smooth static ICM ram pressure is 
proportional to the ICM density picM and the square of 
the galaxy velocity with respect to the Virgo cluster fgai. 
If the ICM is moving with respect to the cluster mean 
velocity the expression for ram pressure yields: 

Pram = PlGuiVgaX - 'UICm)^ , (2) 

where vicm is the galaxy's velocity vector with respect 
to the Virgo cluster, and Wgai is the galaxy velocity with 
respect to the Virgo cluster. 

As noted in Kenney et al. (2004) the orbit of NGC 4522 
is not easy to understand. Its high radial velocity (~ 



1300 kms^^) with respect to the cluster mean velocity 
excludes a classical radial orbit. At its projected distance 
from the cluster center (3.3° ~ 1 Mpc) the density of the 
intracluster medium is an order of magnitude too low to 
strip the galaxy's ISM up to a galactic radius of 3 kpc if 
one assumes a velocity of the galaxy of 1500 km within 
the cluster (Kenney et al. 2004). Kenney et al. (2004) and 
Vollmer et al. (2004) propose two possible scenarios en- 
hancing the ram pressure by a factor of 10: (i) the ICM 
is moving opposite to the trajectory of NGC 4522 due to 
the infall of the M49 group into the Virgo cluster and (ii) 
NGC 4522 is not bound to the cluster and crosses it with 
a velocity of ~ 4000 kms^^ only once. In both cases we 
expect a strongly peaked temporal ram pressure profile 
similar to that for radial orbits: 

Pram — i^max i2 i /2 ' V z 

where tnw is the width of the profile (Vollmer et al. 
2001). Since the observations of Kenney et al. (2004) and 
Vollmer et al. (2004) suggest that ram pressure is close 
to maximum, we only consider model snapshots close 
to peak ram pressure Pmax, which has to be at least 
1000 cm~'^(kms~^)^ to strip gas with a column density 
of 10 M0pc~^ at the observed stripping radius of 3 kpc 
(Kenney et al. 2004). We do not consider higher ram 
pressure maxima, because this value is already a factor 
10 above the estimate based on a static smooth intra- 
cluster medium and a velocity of ~ 1000 kms~^. The 
Gaussian half width tnw is not critical as long as it is 
not too large. Since the velocity of NGC 4522 with re- 
spect to the intracluster medium (flowing or not) is high, 
a value between 50 and 100 Myr is reasonable. Thus, we set 
p,„^^=2000 cm-3(kms-i)2 and iHW=80 Myr. We define 
t = Myr as the time when ram pressure is maximum. 

The ram pressure efficiency also depends on the incli- 
nation angle i between the galactic disk and the ICM wind 
direction (Vollmer et al. 2001). Since the Hi observations 
indicate that stripping occurs more face-on, we do not 
consider edge-on stripping. We made 3 simulations with 3 
different inclination angles between the galaxy's disk and 
the ICM wind direction: (i) i = 45°, (ii) i = 60°, and (hi) 
i = 75°. An inclination of i = 0° means that the galactic 
disk is parallel to the ICM wind direction. 

The last free parameter, the azimuthal viewing angle, 
is chosen in a way to fit the observed Hi distribution and to 
reproduce the positive line-of-sight component of the wind 
direction (the galaxy is moving away from the observer). 
In Figs, m we show the components of the 3D ICM wind 
direction (left panels) and the projected wind direction as 
a function of the azimuthal viewing angle a (right pan- 
els). The maximum of the line-of sight component of the 
ICM wind is largest for i = 45° and smallest for i — 75°. 
The wind components along the right ascension and dec- 
lination are all negative, i.e. the ICM wind blows from 
the south east. Whereas the two components are roughly 
equal for i = 45° and i = 60°, the RA components is 
smaller than the DEC components for i — 75°. As the 
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Fig. 1. Inclination angle between the galactic disk and the 
ICM wind i = 45°, i = 60°, and i = 75°. Left panels: com- 
ponents of the 3D ICM wind direction as a function of the 
azimuthal viewing angle a. Right panels: projected wind 
direction with respect to the galaxy's minor axis (counted 
counter-clockwise) as a function of a. The arrows indicate 
the best fit model a — 350° and the azimuthal viewing 
angle for a projected wind angle of —15°. 



maximum of the line-of sight component the maximum of 
the projected wind direction also varies with i. It is max- 
imum for i — 45° and minimum for i = 75°. In all cases, 
an azimuthal viewing angle of a = 350° leads to a positive 
projected wind direction smaller than 6°. However, if one 
assumes a given projected wind direction of —15°, the cor- 
responding azimuthal viewing angles vary: a (45°) = 20°, 
a(60°) = 30°, a(75°) = 90°. 

None of the simulations shows significant extraplanar 
gas before t = Myr. In Fig.[21we show 3 timesteps Myr, 
50 Myr, and 100 Myr after the ram pressure maximum. 
The simulation with an inclination angle between the 
galaxy's disk and the ICM wind direction of i = 45° does 
not lead to the observed stripping radius neither shows 
significant extraplanar gas. The simulation with i = 75° 
(Fig.[5^-c) leads to the correct stripping radius and shows 
extraplanar gas at t = 50 Myr (Fig. [2h) and t = 100 Myr 
(Fig. 13). However, the column density of this extraplanar 
gas is smaller compared to the Hi observations. In partic- 
ular, the observed detached south western extraplanar Hi 
maximum is missing. Only the simulation with i — 60° 
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Fig. 3. Snapshot at i = 40 Myr of the simulation with 
i = 60°. The azimuthal viewing angle is a = 30°. Left 
panel: gas distribution. Right panel: velocity field. 



at t = 50 Myr (Fig. ^) can reproduce the observed Hi 
distribution. About 50 Myr later the column density of 
the northern extraplanar gas has already decreased signif- 
icantly and the shape of the southern extraplanar gas has 
changed (Fig.gf). 

As a conclusion, we rule out the i = 75° simulations, 
because the galaxy's small line-of-sight velocity compo- 
nent (Fig. ^ would lead to galaxy and/or ICM velocities 
which are beyond any reasonable values (> 4000 kms^^). 
We also rule out the i — 45° simulations, because it can 
not reproduce the observed extraplanar gas. We are thus 
left with the i — 60° simulation. Recent HST observations 
of NGC 4522 (Kenney et al., in prep.) yield a projected 
wind direction of ~ —15°. The corresponding azimuthal 
viewing angle is a = 30°. Using this a we realized that 
the observations can be best reproduced with t — 40 Myr 
instead of i = 50 Myr used before (Fig. (SJ. Whereas the 
Hi gas distribution is reproduced in a satisfactory way, the 
velocity isocontours of the western extraplanar gas are dif- 
ferent from the observed ones (Fig. IHJ. 

As an additional test we approximated the temporal 
ram pressure profile with a step function, i.e. the galaxy is 
moving abruptly into a region of high intracluster medium 
density and/or velocity. The step occurs a,t t — Myr and 
the subsequent ram pressure is assumed to be constant at 
a level of the maximum ram pressure of the Lorentzian 
(Eq.OJ. Fig. H shows the ISM evolution for i = 60°. None 
of the timesteps show the observed distribution of the ex- 
traplanar high column density gas (Fig. |Hl . We therefore 
exclude a constant ongoing ram pressure. 

Therefore, our best fitting model has the following 
parameters (see Eq. OJ: Pmax = 2000 cm~'^(kms~^)^, 
iHW=80 Myr, i = 60°, and t = 50 Myr. The associated 
uncertainties are Ai ~ 10°, Aa ~ 20°, and At ^ 10 Myr. 
Subsequently, the radial velocity of NGC 4522 represents 
33 % of its total velocity within a static ICM. The total 
velocity of NGC 4522 with respect to a static ICM is thus 
'yN4522 ^ 3500 kms~^. Together with a local ICM density 
of nicM = 10~^ kms~^ this is consistent with our assumed 
peak ram pressure of 2000 cm~^(kms~^)^ within a factor 
of 1.6. 
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Fig. 2. Model snapshots for different simulations. The ram pressure profile for all simulations is given by Eq. The 
inclination angle i between the orbital and the galaxy's disk plane and the timestep of the snapshot are varied. Left 
column: t — Myr. Middle column: t = 50 Myr. Right column: t — 100 Myr. Upper row: i = 45°. Middle row: i = 60°. 
Lower row: i = 75° . 



The large scale evolution of this NGC 4522 simula- 
tion can be seen in Fig. [3 The ICM wind begins to drive 

out the ISM of NGC 4522 only at 100 Myr. The 

gas is expelled to the west. At the first stage of stripping 
{t ~ —30 Myr) the extraplanar gas of highest surface den- 
sity is found in the south-west. At < ~ 10 Myr both, the 
south-western and north-western parts of the extraplanar 
gas, show a relatively high column density. After the ram 
pressure maximum has occurred (t = Myr) the surface 
density of the extraplanar gas west to the galaxy center 
has significantly increased. 



According to these simulations we observe the galaxy 
~ 50 Myr after the ram pressure maximum. In the two 
different scenarios described above this means that (i) 
the galaxy encountered the intracluster medium of high- 
est density ~ 50 Myr ago or (ii) the galaxy just leaves 
the region of maximum ICM velocity. The final gas dis- 
tribution is shown in Fig. IHl The extraplanar gas has a 
complex asymmetric three-dimensional structure. A part 
of the ISM initially located at galactic radii between 3 and 
5 kpc is stripped in a ring-like configuration, located up 
to ~ 4 kpc above the stellar disk. 
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Fig. 6. Comparison snapsliot of the gas distribution of NGC 4522. Left panel: view of the disk using the observed 
position angle and inclination of NGC 4522. Right panel: the disk seen face-on. The arrow indicates the direction of 
the wind, i.e. opposite to the galaxy's motion. 




Fig. 4. Snapshots of the simulation using a step function 
for the temporal ram pressure profile with i = 60°. The 
azimuthal viewing angle is a = 350°. 



4. Comparison with VLA Hi observations 

In this section the model snapshot of Fig. is compared 
with the VLA Hi observations of Kenney et al. (2004). 
For this purpose we assume that only gas with a volume 
density greater than ~ 10 cm^'^ is in form of neutral hy- 
drogen. This corresponds to the high end of the densities 
of the warm ionized medium (Boulares & Cox 1990). The 
lower density gas is assumed to be ionized. The compari- 
son between the model and observations is done based on 
moment maps (Sect. 14. 1|) and a position- velocity diagram 
(Sect. EH). 



4.1. The moment maps 

Since ram pressure selectively affects the gas, the model 
and observed stellar disks are symmetric whereas the 
model and observed gas distributions are highly asym- 
metric (Fig. Cj), i.e. the ISM is pushed to the west of the 
galactic disk. 

The Hi distribution (moment 0), velocity field (mo- 
ment 1), and velocity dispersion (moment 2) are shown in 
Fig. |S1 The following observed properties of the Hi data 
are reproduced: 

— the model Hi deficiency is 0.5 compared to the ob- 
served value of 0.6 ± 0.2, 

— the gas disk is truncated at a radius of about 3 kpc, 

— within the disk more gas is found in the northern part 
than in the southern part, 

— the existence of a high surface density extraplanar gas 
component to the west of the galactic disk, 

— the extraplanar gas distribution shows two maxima; 
the south-western maximum is more prominent than 
the north-western one, 

— the velocity field of the extraplanar gas is relatively 
regular, 

— the extraplanar gas has less extreme velocities than the 
nearby disk emission, i.e. the velocity contours curve 
away curve away from the minor axis. 

Thus, the overall agreement between the model and ob- 
served first two moments (gas distribution and velocity 
field) is satisfactory. Whereas 40% of the observed Hi mass 
resides in the extraplanar component, the model yields a 
fraction of 65%. This difference is likely due to the un- 
known difference between the initial gas distribution of 
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Fig. 5. Evolution of the model stellar (1st column) and 
gas disk (2nd column). The major axis position angle and 
inclination of NGC 4522 are PA = 214° and « = 80°, re- 
spectively. The arrow indicates the direction of ram pres- 
sure, i.e. it is opposite to the galaxy's velocity vector, and 
its size is proportional to pWgai- Maximum ram pressure 
occurs at t = Myr. The timestep of each snapshot is 
marked in each panel showing the stellar disk. 



the model and/or due to the fact that we detect only a 
fraction of the stripped gas due to ISM expansion, heating 
and ionization. 

The model and observed velocity dispersions are dif- 
ferent. In the observations, the peak Hi dispersion is not 
at the nucleus, but in the extraplanar gas, ~ 30" west of 
the nucleus. Even within the disk, the maximum in the 
observed dispersion is offset from the nucleus by ~ 10". 
Our model shows a maximum in the center, a region of 




Fig. 7. Stellar (contours) and gas (greyscale) distribution 
of NGC 4522. Left panel: model. Right panel: Hi observa- 
tions (Kenney et al. 2004). 

marginally larger linewidth to the north, and a region 
of larger linewidth within the south western extraplanar 
gas. The central model peak is due to the model rotation 
curve which is somewhat steeper than the observed rota- 
tion curve. Thus, although we find a small enhancement 
of the linewidth in the south western extraplanar gas, its 
linewidth is a factor of 3 smaller than the observed one. 



4.2. Position-velocity diagram 

The position-velocity diagram along the major axis for the 
disk and extraplanar regions are shown in Fig.|5| The disk 
region is defined as the region between —20" and -1-11" 
from the major axis. The extraplanar region is defined as 
the region between -1-11" and -1-50" from the major axis. 
The following observed characteristics are reproduced by 
the model: 

— the disk gas kinematics are relatively symmetric and 
normal, despite the strong north-south asymmetry of 
the gas distribution, 

— the extraplanar gas has lower (blueshifted) velocities 
with respect to the disk gas. 

On the other hand, our model does not reproduce the 
observed large linewidth (up to 150 kms^^ FWZI). 



5. MHD simulations 

5.1. The model 

Otmianowska-Mazur & VoUmer (2003) studied the evolu- 
tion of the large scale magnetic field during a ram pres- 
sure stripping event. They calculated the magnetic field 
structure by solving the induction equation on the veloc- 
ity fields produced by the dynamical model. The polarized 
radio continuum emission has been calculated by assum- 
ing a Gaussian spatial distribution of relativistic electrons. 
This procedure allowed them to study the evolution of the 
observable polarized radio continuum emission during a 
ram pressure stripping event. 
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We apply the same procedure as Otmianowska-Mazur 
& VoUmer (2003) on a similar ram pressure stripping event 
(Sect. 13). The Zeus3D code (Stone & Norman 1992a and 
b) is used to solve the induction equation: 

dB / dt ^ mi{v X B) - T:oi{ri TotB) (4) 

where B is the magnetic induction, v is the large-scale ve- 
locity of the gas, and r] is the coefHcient of a turbulent dif- 
fusion. We use a physical diffusion of = 5 x 10^^ cm^s^^ 
(Elstner et al. 2000). The estimated numerical diffusion 
is an order of magnitude smaller and thus does not affect 
our simulations. We do not use an a dynamo. The initial 
magnetic field is purely toroidal with a strength of 10 (iG. 

The induction equation is solved on a rectangular co- 
ordinates {XYZ). The number of grid points used is 
171x171x71 along the X, Y and Z axis, respectively. This 
corresponds to the grid spacing of 200 pc in the galactic 
plane and of 300 pc in the Z direction, resulting in a size 
of the modeled box of 34.2 kpc x 34.2 kpc x 21.3 kpc. 
Since the N-body code is discrete whereas the MHD code 
is using a grid, we have to interpolate the discrete ve- 
locities on the grid. This is done using a method known 
as "Kriging" with a density-dependent smoothing length 
(Isaaks & Srivastava 1989). It turned out that we had to 
use a large smoothing length to suppress the noise in the 
velocity field of the outer disk, which is due to a small, 
local particle density. In this way we avoid numerical ar- 
tifacts of the magnetic field distribution at the outer disk. 
As the spline interpolation used in Otmianowska-Mazur & 
Vollmer (2003), this has the consequence that the velocity 
field at the edge of the gas distribution is more extended 
than the gas distribution itself. Since there are gradients 
in this velocity field, induction leads to a magnetic field 
which extends beyond the edge of the gas distribution. 
This affects the polarized emission beyond the gas distri- 
bution, but not inside, which is what we are interested 
in. 

The evolution of the polarized radio continuum emis- 
sion without Faraday rotation is presented in Fig. ^| The 
timesteps are the same as in Fig.|Sl The disk rotates clock- 
wise. The gas surface density, which is smoothed to a res- 
olution of ^ 100 pc, is shown in greyscales, the polarized 
radio continuum emission as contours, and the magnetic 
field vectors projected on the plane of the sky as lines. We 
assume a Gaussian distribution of relativistic electrons in 
R and z directions: n^ci = n-o exp ( — (r/rpy^)^), where 
rpi, = 5 kpc and — 0.5 kpc. This translates into a 
FWHM of 8.3 kpc and 1.2 kpc, respectively. By assuming 
this smooth distribution we imply no equipartition be- 
tween the energy densities of total cosmic rays and total 
magnetic field. The galaxy moves to the lower left corner, 
i.e. the ram pressure wind comes from the this direction. 
The ram pressure maximum occurs a,t t — Myr. As al- 
ready seen in Sect. El ram pressure begins to push the gas 
to the upper right corner at t ^ —100 Myr. The straight 
cut-off of the gas distribution in the upper right quadrant 
is due to the edge of our computational volume. 



At t > —100 Myr the gas is compressed at the lower 
left side of the galaxy. At the same time most of the gas 
initially located in the outer disk (i? > 3 kpc) is pushed 
out of the disk plane (Fig. (S)). At t > Myr a promi- 
nent maximum of the polarized radio continuum emis- 
sion forms in the lower left quadrant. Due to our large 
density-dependent smoothing length, a part of the large- 
scale magnetic field in the lower left quadrant of the galaxy 
does not follow the compression and stays at large galac- 
tic radii in the disk plane (see also Otmianowska-Mazur 
& Vollmer 2003). Therefore, the distribution of polarized 
intensity extends beyond that of the gas which is mostly 
due to our smoothing algorithm. However, a decoupling 
of the large-scale magnetic field from the gas flow is not 
entirely excluded (see .e.g. Soida et al. 2001). Recently, 
Beck et al. (2005) could show that the regular magnetic 
field can decouple from the cold gas. The maximum of 
polarized emission within the gas disk in the lower left 
direction is real and due to gas compression. As observed 
in Otmianowska-Mazur & Vollmer (2003) after maximum 
compression at t = Myr this inner maximum is taken 
along with rotation towards the upper left part of the im- 
age {t > 20 Myr). At t = 50 Myr it is located to the left 
of the galaxy center, at t = 90 Myr it is located on the 
upper left part of the gas disk. Since the FWHM of the 
relativistic electron distribution in the vertical direction is 
1.2 kpc, the magnetic field taken away vertically with the 
gas is no longer observable in polarized radio continuum 
emission because of the low density of relativistic electrons 
there. 

The evolution of the polarized radio continuum 
emission (Fig. I10|) is qualitatively similar to that of 
Otmianowska-Mazur & Vollmer (2003), where an inclina- 
tion angle between the disk and the ICM wind direction 
of « = 20° was assumed. Thus, the Otmianowska-Mazur & 
Vollmer (2003) simulation corresponded to a more edge-on 
stripping, whereas in our simulation ram pressure strip- 
ping occurs more face-on {i = 60°). Consequently, the 
gas compression is less pronounced in the present case. 
However, as in the Otmianowska-Mazur & Vollmer (2003) 
simulations we observe that shortly after the timestep of 
maximum ram pressure the maximum of polarized ra- 
dio continuum emission within the compressed gas disk 
is moving to the upper part of the image due to galac- 
tic rotation. However, ram pressure is still quite high and 
the ISM is still leaving the galactic disk at t = 50 Myr. 
This is preferentially happening on the side where galactic 
rotation and the ICM wind are parallel, i.e. in the north 
east. 

As a conclusion, the polarized radio continuum emis- 
sion outside the gas disk at timesteps t > —20 Myr is most 
probably due to our large density-dependent smoothing 
length due to the intrinsic difficulties of interpolating and 
extrapolating the 3D velocity field over a sharp edge of 
the gas distribution where no particles are found beyond 
a certain radius. In this case an interpolation and extrap- 
olation with a density dependent smoothing length can al- 
ways give rise to artifical velocity gradients which lead to 
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an enhancement of the magnetic field. The here employed 
"Kriging" method with a density dependend smoothing 
length turned out to be a good compromise between the 
need of a smooth velocity and a minimum of artificial ve- 
locity gradients at the edge of the gas distribution. 

However, the maximum of polarized radio continuum 
emission inside the compressed gas disk is real. Since we 
know the cause of polarized emission extending beyond the 
gas disk and since this feature and the real inner maximum 
of polarized emission are very close in space and since we 
did not want to modify the model by hand with a risk of 
introducing uncontrollable artifacts, we did not attempt to 
remove the part of the large-scale magnetic field extending 
beyond the gas disk. We interpret the observed maximum 
of polarized radio continuum emission within the gas disk 
as being due to gas compression. 

5.2. Comparison with VLA polarized radio continuum 
emission 

The distribution of polarized radio continuum emission is 
now projected using the position and inclination angles of 
NGC 4522, t — 50 Myr, and the azimuthal angle as in 
Sect.m 

The resulting model polarized radio continuum emis- 
sion distribution and the model magnetic field vectors on 
the model gas distribution are shown in Fig. 1111 together 
with our 6 cm VLA observations (VoUmer et al. 2004). 
As our observations, the model shows a polarized radio 
continuum maximum to the east of the galaxy center and 
a weaker extension to the north. The observed offset of 
the polarized emission maximum to the east translates 
the fact that the inclination angle of the galaxy is such 
that the eastern side of the disk is the compressed side 
which is the far side of the galaxy. Since the polarized 
emission extending beyond the edge of the gas distribu- 
tion and the real compressed large-scale magnetic field 
are contiguous, the shown distribution of polarized radio 
continuum emission reflects that of the compressed com- 
ponent. The observed elongation of the polarized radio 
continuum emission can be explained as follows: due to 
gas compression a polarized emission maximum is created 
in the south-cast. Shortly after the occurrence of the ram 
pressure maximum, the large-scale magnetic field follows 
rotation to the north-east. Thus, the observed polarized 
radio continuum emission corroborates our previous find- 
ing (based on the comparison between the model and the 
Hi data) that we observe NGC 4522 after the occurrence 
of maximum ram pressure (Sect.|21). 

6. Discussion 

6.1. Stripping of the cold and warm Hi 

For face-on stripping, our simulations are consistent with 
the Gunn & Gott (1972) ram pressure estimate. Our dy- 
namical model reproduces qualitatively and quantitatively 
the observed Hi gas distribution, velocity field (Fig.lHjl and 



polarized radio continuum emission distribution fFig.lll|). 
However, it fails in reproducing the observed large Hi 
linewidths in the extraplanar component (Fig.lH)). This is 
most probably due to a lack of gas physics in the model. 
The large observed linewidth of ~ 100 kms~^ cannot be 
thermal, because this would translate into temperatures 
where the gas is ionized. Thus, they are due to directed, 
chaotic or turbulent gas motions. Since the Hi line profiles 
tend to smaller velocities, i.e. to velocities closer to that of 
the Virgo cluster, one possible explanation is that a part 
of the atomic gas is stripped more efficiently than the rest. 
Since the column density of the low velocity gas is small, 
we speculate that it is a low density warm (T ~ 8000 K) 
atomic gas phase which is stripped more efficiently. In this 
scenario the large linewidth is due to directed motions. 
In a second scenario ram pressure leads to shocks in the 
interstellar medium which increase their chaotic or turbu- 
lent velocity dispersion dramatically. Again, as the column 
density of the low velocity gas is small, this might only 
concern the low density warm component of the atomic 
gas which has a larger volume filling factor than the cold 
atomic gas (T ~ 100 K). Both scenarios can explain the 
observed large Hi linewidths in the extraplanar compo- 
nent. In both scenarios it is the warm diffuse atomic gas 
phase which produces the large linewidths. 

6.2. The stripping efficiency 

As already mentioned in Sect. |21 our knowledge about the 
orbit of NGC 4522 within the Virgo intracluster medium 
is poor. Its high radial velocity with respect to the clus- 
ter mean and its large distance from the cluster center 
exclude a simple radial orbit. Kenncy et al. (2004) esti- 
mated ram pressure due to a static intracluster medium 
at the projected location of NGC 4522. Applying the Gunn 
& Gott (1972) criterion and assuming an ICM density of 
"-ICM = 10"'' cm~'^ and a galaxy velocity of 1500 kms~^, 
they found that the gravitational restoring force for a gas 
cloud with a surface density of 10 Mopc~^ is an order 
of magnitude higher than the force due to ram pressure. 
This led them to the conclusion that ram pressure is higher 
than the "standard value" . This could in principle be due 
to a higher galaxy velocity or a higher stripping efficiency ^ 
withpiam — £.p{vga\^vicM.)^ (sec e.g. Roediger & Hensler 
2005). 

Our simulation yields a projected ICM wind direc- 
tion to the north east which is roughly aligned with the 
galaxy's minor axis (6° ± 20°; Fig. EJ. This is corrobo- 
rated by the recent results of Kenney et al. (in prep.) 
analysing dust extinction features in deep HST images 
of NGC 4522. They found gas clouds beyond the strip- 
ping radius in the direction of the ICM wind. Beyond 
the main gas/dust truncation radius, they found several 
elongated dust clouds, all with similar position angles of 
— 17° ±3° with respect to the minor axis, likely indicating 
the projected ICM wind direction. Moreover, our simula- 
tions yield an inclination angle between the disk and the 
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ICM wind direction of i = 60° ± 10° (more face-on strip- 
ping). Thus the wind angle is relatively face-on despite 
the fact that the galaxy us nearly edge-on and has a high 
line-of-sight velocity, implying that there must be a large 
component of the ICM wind in the plane of the sky. Based 
on our model the total velocity of the galaxy relative to a 
static ICM is - 3500 kms^^ (see Sect.E)). With this ve- 
locity the Gunn & Gott criterion is approximately fulfilled 
and the stripping efficiency is close to one. 

6.3. Static or moving ICM 

The comparison between Hi observations and our model 
(Sect.^ suggests that ram pressure acting on NGC 4522 
is close to its maximum but has already passed it. The 
comparison between the observed and modeled polarized 
radio continuum emission confirms that the maximum has 
occurred recently (~ 50 Myr ago), because the polarized 
radio continuum emission extends to the north-east, in the 
direction of galactic rotation. 

Thus, if the galaxy is crossing the cluster with a very 
high velocity (~ 4000 km^^) and the ICM is static and 
smooth, NGC 4522 has just passed the point of highest 
intracluster medium density, i.e. it has just passed the 
point of smallest distance to M87. On the other hand, if 
the intracluster medium is moving due to the infall of the 
M49 group of galaxies as suggested by Kenney et al. (2004) 
and VoUmer et al. (2004), the galaxy has just passed the 
region of highest intracluster medium velocities. In this 
case ram pressure is given by Eq.[21 If the main component 
of fgai is radial and thus Wgai ~ 1500 kms~^, the ICM 
velocity vector vicm must have dominant components in 
the plane of the sky. In the following we give an example: 

We define the signs of velocity vectors as positive to 
the west, north, and away from the observer. Assuming a 
galaxy velocity vector of v^^i =(-675, -675, 1155) kms~^ 
and taking the total velocity vector from our simula- 
tion i^tot ~ 'Ugai — i'iCM=(-2560, -2090, 1155) kms~^ 
leads to a galaxy velocity of Ugai = 1500 kms~^ and 
vtot = 3500 kms~^ (see Sect.EJ. The resulting intraclus- 
ter medium velocity is then t^icM = fgai — 'i'tot = (1885, 
1415, 0) kms~^ and the total ICM velocity is vicm = 
2360 kms~^. Thus the motion of the intracluster medium 
which might be partly due to the infall of the M49 sub- 
group is towards the north west and exclusively in the 
plane of the sky. This ICM velocity is 1000 kms~^ higher 
than that derived from X-ray data assuming pressure equi- 
librium (Shibata et al. 2001). Alternatively, a local ICM 
density enhancement of a factor ~4 together with an ICM 
velocity of ~1000 kms~^ can account for the observed 
stripping radius of NGC 4522. This is consistent with the 
radial velocity of M49 (~ 1000 kms~^) which is very close 
to the velocity of the Virgo cluster. The derived direction 
of this putative ICM flow might also be related to the X- 
ray detection of a gas compression region in the north of 
M49 (Irwin & Sarazin 1996, Biller et al. 2005) 



7. Conclusions 

We present a dynamical model for the evolution of the 
interstellar medium of the Virgo cluster spiral galaxy 
NGC 4522, which presently undergoes strong ram pressure 
stripping (Kenney et al. 2004). We confront the observed 
with the simulated moment maps. In addition, we solve 
the induction equation on the velocity fields of the dy- 
namical model to calculate the large-scale magnetic field. 
Assuming a smooth relativistic electron distribution we 
obtain the model distribution of polarized radio contin- 
uum emission which can be directly compared to observa- 
tions (VoUmer et al. 2004). 
We conclude that 

— the model successfully reproduces the observed Hi gas 
distribution and velocity field; 

— the model fails in reproducing the observed large Hi 
linewidths (^ 100 kms^^); 

— the MHD model successfully reproduces the observed 
polarized radio continuum emission; 

We suggest that large observed Hi linewidths is due to 
the warm diffuse component of the atomic gas which is ei- 
ther stripped more efficiently or whose velocity dispersion 
is increased by ram pressure induced shocks. 

The model confirms that NGC 4522 undergoes ram 
pressure which is close to its maximum. This maximum 
ram pressure of ^ 2000 cm^'^(kms~^)^ occurred about 
50 Myr ago. The inclination angle between the galaxy's 
disk and the ICM wind direction is i = 60° ± 10° and the 
ICM wind blows from the south-east. The projected wind 
angle with respect to the galaxy's minor axis is 6° ± 20°. 
Since the galaxy is located at a projected distance of 
~ 1 Mpc from the cluster center where the intracluster 
medium density is by far too low to cause significant strip- 
ping if it is static, NGC 4522 is either on an unbound orbit 
within the cluster and just passed the region of highest 
Virgo intracluster medium density or, if the intracluster 
medium is moving, the galaxy just leaves the region of 
highest intracluster medium velocities. If the intracluster 
medium is static and smooth the galaxy has a total veloc- 
ity of about 3500 kms^^ with respect to the Virgo cluster. 
If the intracluster medium is moving possibly due to the 
infall of the M49 group of galaxies, we derive a velocity 
of the infalling intracluster medium of about 2400 kms~^ 
with respect to the Virgo cluster. In this case the intra- 
cluster medium associated with the M49 group is moving 
towards the north east with a negligible radial velocity 
component. 

This study shows the strength of combining high res- 
olution Hi and polarized radio continuum emission with 
detailed numerical modeling of the evolution of the gas 
and the large-scale magnetic field. 
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Fig. 8. Moment maps of NGC 4522. Left panels: model. Right panels: Hi observations (Kenney et al. 2004). The same 
contouring scheme is used for both model and data. Upper panels: gas distribution; contours spaced by \f2. Middle 
panel: velocity field; contours spaced 10 kms"-^. Lower panel: velocity dispersion; contours: 15, 20, 25, 30, 35, 40, 
45 kms~^. 
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Fig. 9. Major axis position-velocity diagram for the disk (greyscale) and extraplanar (contours) component. Left panel: 
model. Right panel: Hi observations (Kenney et al. 2004). Contour levels are 1, 2, 3, 4, 5, 6 times the peak value. 
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Fig. 10. Evolution of the polarized radio continuum emission. Time steps are the same as in Fig. |S1 The disk is 
seen face-on and rotates clockwise. Contours: polarized radio continumn emission. The magnetic field vectors are 
superimposed onto the gas surface density. The timestep is indicated on top of each snapshot. The assumed beamsize 
for the polarized radio continumn emission is 20". The arrow in the last panel indicates the direction of the ram 
pressure wind. 
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Fig. 11. Greyscale: Hi gas distribution with a resolution of -^20". Contours: polarized radio continuum emission. Left 

panel: model. The FWHM of the rclativistic electron distribution in the radial and vertical directions are 8.3 kpc and 
1.2 kpc, respectively. Right panel: observations (VoUmer et al. 2004). 



